A two way fully coupled fluid structure simulation of human ureter peristalsis.

Numerical simulations of ureter peristalsis have been carried out in the past to understand both the flow field and ureter wall mechanics. The main objective of the current investigations is to have a better understanding of the urine transport due to the peristalsis in the ureter, thus making the information helpful for a better treatment and diagnosis of ureteral complications like urine reflux. In the current study, a numerical simulation is performed using a finite-element-based solver with a two-way fully coupled fluid structure interaction approach between the ureter wall and urine. For the first time, the ureter wall is modeled as an anisotropic hyper-elastic material based on experiments performed in previous literature on the human ureter. Peristalsis in the ureter is modeled as a series of isolated boluses. By observing the flow field it is clear that the peristalsis mechanism has a natural tendency to create a backflow as the isolated bolus moves forward. As a result, the urine can flow back from the bladder to the ureter at the ureterovesical (ureter-bladder) junctions, if the one-way valve starts to malfunction.

A three-dimensional (3D) two-way coupled fluid-structure interaction (FSI) study of peristaltic flow in obstructed ureters.

Obstruction in the ureter flow path is one of the most common problems in urinary-related diseases. As the ureter transports the urine using the expansion bolus created by the peristaltic pulses, an obstruction in its path can cause unwanted backflow and can also result in damage to the wall. But in order to understand this further, and specifically to quantify and parametrize the effect of the obstruction in the ureter, a detailed study investigating various level of obstructions in peristaltic ureter flow is necessary. In the current study, full 3D numerical simulations of peristalsis in an obstructed ureter are carried out using a finite element solver along with a two-way coupling between the fluid and structural domain with the arbitrary Eulerian-Lagrangian method. Analysis of the results shows that the larger the obstruction, the higher the wall shear stress and pressure gradient in the fluid. In addition, the amount of backflow increases with increase in the obstruction.

A fluid-structure interaction (FSI)-based numerical investigation of peristalsis in an obstructed human ureter.

Urine moves from the kidney to the bladder through the ureter. A series of compression waves facilitates this transport. Due to the highly concentrated mineral deposits in urine, stones are formed in the kidney and travel down through the urinary tract. While passing, a larger stone can get stuck and cause severe damage to ureter wall. Also, stones in the ureter obstructing the urine flow can cause pain and backflow of urine which in turn might require surgical intervention. The current study develops a 2D axisymmetric numerical model to gain an understanding of the ureter obstruction and its effects on the flow, which are critical in assessing the different treatment options. Transient computational analysis involving a two-way fully coupled fluid-structure interaction with the arbitrary Lagrangian-Eulerian method between the ureteral wall and urine flow is conducted with an obstruction in the ureter. The ureter wall is modeled as an anisotropic hyperelastic material, data of which, is based on biaxial tests on human ureter from previous literature, while the incompressible Navier-Stokes equations are solved to calculate urine flow. A finite element-based monolithic solver is used for the simulations here. The obstruction is placed in the fluid domain as a circular stone at the proximal part of the ureter. One of the objectives of this study is to quantify the effect of the ureteral obstruction. A sharp jump in pressure gradient and wall shear stress, as well as retrograde urine flow, is observed as a result of the obstruction.

Three-Dimensional Numerical Simulations of Peristaltic Contractions in Obstructed Ureter Flows.

Ureteral peristalsis can be considered as a series of waves on the ureteral wall, which transfers the urine along the ureter toward the bladder. The stones that form in the kidney and migrate to the ureter can create a substantial health problem due to the pain caused by interaction of the ureteral walls and stones during the peristaltic motion. Three-dimensional (3D) computational fluid dynamics (CFD) simulations were carried out using the commercial code ansys fluent to solve for the peristaltic movement of the ureter, with and without stones. The effect of stone size was considered through the investigation of varying obstructions of 5%, 15%, and 35% for fixed spherical stone shape. Also, an understanding of the effect of stone shape was obtained through separate CFD calculations of the peristaltic ureter with three different types of stones, a sphere, a cube, and a star, all at a fixed obstruction percentage of 15%. Velocity vectors, mass flow rates, pressure gradients, and wall shear stresses were analyzed along one bolus of urine during peristalsis of the ureteral wall to study the various effects. It was found that the increase in obstruction increased the backflow, pressure gradients, and wall shear stresses proximal to the stone. On the other hand, with regard to the stone shape study, while the cube-shaped stones resulted in the largest backflow, the star-shaped stone showed highest pressure gradient magnitudes. Interestingly, the change in stone shape did not have a significant effect on the wall shear stress at the obstruction level studied here.

Magnetohydrodynamic simulations of hypersonic flow over a cylinder using axial- and transverse-oriented magnetic dipoles.

Numerical simulations of magnetohydrodynamic (MHD) hypersonic flow over a cylinder are presented for axial- and transverse-oriented dipoles with different strengths. ANSYS CFX is used to carry out calculations for steady, laminar flows at a Mach number of 6.1, with a model for electrical conductivity as a function of temperature and pressure. The low magnetic Reynolds number (<<1) calculated based on the velocity and length scales in this problem justifies the quasistatic approximation, which assumes negligible effect of velocity on magnetic fields. Therefore, the governing equations employed in the simulations are the compressible Navier-Stokes and the energy equations with MHD-related source terms such as Lorentz force and Joule dissipation. The results demonstrate the ability of the magnetic field to affect the flowfield around the cylinder, which results in an increase in shock stand-off distance and reduction in overall temperature. Also, it is observed that there is a noticeable decrease in drag with the addition of the magnetic field.

Nonoscillatory central schemes for hyperbolic systems of conservation laws in three-space dimensions.

We extend a family of high-resolution, semidiscrete central schemes for hyperbolic systems of conservation laws to three-space dimensions. Details of the schemes, their implementation, and properties are presented together with results from several prototypical applications of hyperbolic conservation laws including a nonlinear scalar equation, the Euler equations of gas dynamics, and the ideal magnetohydrodynamic equations. Parallel scaling analysis and grid-independent results including contours and isosurfaces of density and velocity and magnetic field vectors are shown in this study, confirming the ability of these types of solvers to approximate the solutions of hyperbolic equations efficiently and accurately.